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SUMMARY 


The  research  on  boiling  in  confined  spaces  has  been  completed. 

This  problem  is  of  great  importance  to  the  boiling  induced  corrosion  in 
the  steam  generator  crevices  between  the  tube  and  the  support  structures. 

In  the  report  of  1981,  analysis  of  single  phase  flow,  two  phase  flow,  and 
dryout  in  crevices  have  been  presented.  Experimental  results  of  boiling 
and  dryout  in  vertical  crevices  with  closed  bottom  were  also  reported. 

In  the  report  of  1982,  the  results  of  forced  convective  boiling  heat 
transfer  and  critical  heat  flux  in  crevices  are  shown.  In  this  report, 
the  boiling  heat  transfer  and  critical  heat  flux  in  horizontal  crevices 
are  reported. 

The  research  on  cross  flow  boiling  in  tube  bundles  proceeds  smoothly. 

In  this  report,  the  analysis  of  single  phase  flow  pressure  drop  and  con¬ 
vective  heat  transfer  in  tube  bundles  at  cross  flow  of  very  low  or  very 
high  Reynolds  numbers  are  presented.  The  experimental  setup  has  been 
completed  and  data  are  taken  presently.  > 


CONCLUSION 


Through  the  systematical  studies  of  the  boiling  heat  transfer  in 
confined  spaces  the  following  general  conclusions  can  be  made: 

(1)  The  boiling  heat  transfer  in  confined  spaces  usually  has  a 
higher  heat  transfer  coefficient  (higher  heat  flux  at  a 
specific  wall  superheat)  than  the  conventional  conditions. 
This  is  due  to  the  effective  thin  film  evaporation  in  the 
crevices  under  the  squeezed  pan-cake  shape  bubbles. 

(2)  The  boiling  regimes  can  be  classified  as  isolated-deformed 
bubbles,  coalesced-deformed  bubbles,  and  the  slightly 
deformed  bubbles.  The  existence  of  the  regimes  is  strongly 
dependent  upon  the  Bond  number  of  the  gap  as  smaller, 
close  to,  or  larger  than  unity.  The  Boiling  number  has 
secondary  order  effect  to  the  boiling  regime. 

(3)  The  critical  heat  flux  usually  occurs  at  highly  oscillatory 
conditions.  The  occurrence  of  CHF  in  pools  can  be  analyzed 
according  to  flooding  or  momentum  balance  criterion. 
Theoretical  modeling  and  semi -empirical  correlations  are 
established.  For  convective  conditions,  the  conventional 
correlation  can  be  used  with  a  correction  factor  applied. 

The  correction  factor,  again,  depends  upon  the  Bond  number 
and  L/D  ratio  of  channel.  The  CHF  decreases  with  decreasing 
gap  sizes.  They  occur  due  to  the  dryout  of  the  thin  liouid 


(4)  With  this  improved  understanding  the  designers  are  able 
to  avoid  the  occurrence  of  the  dryout  and  the  induced 
corrosion  in  the  crevices  which  are  in  adjacent  to  heated 
parts  in  boilers.  A  computer  program  has  been  established 
to  predict  the  dryout  pattern  in  any  geometries.  It  can 
be  used  to  examine  the  performance  of  existing  equipment 
as  well  as  to  prove  any  new  design  concepts. 

The  research  on  cross  flow  heat  transfer  in  tube  bundles  indicated 
that  at  low  speed  creeping  flow  the  flow  streams  are  symmetric  with  a 
cross-sectional  plane;  and,  at  very  high  Reynolds  number  single  phase 
convection,  and  possibly  at  the  convective  boiling  conditions  the  turbu 
lence  is  very  intense  and  homogeneous.  Therefore  the  heat  transfer  on 
tubes  can  be  described  in  a  very  simple  manner  that  contains  mainly  the 
solidity  ratio  (the  cross-sectional  ratio  of  tubes  to  heat  exchangers). 
As  a  result,  the  complicated  Figures  and  Tables  in  presenting  data  of 
tube  bundles  with  various  configuration  can  be  greatly  simplified  with 
short  formulations.  Further  study  is  in  progress  to  confirm  this  under 
standing  on  a  much  broader  data  base. 
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Chapter  1 


Pool  Boiling  Heat  Transfer 
in  Narrow  Horizontal  Annular  Crevices 
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ABSTRACT 


Experimental  results  of  the  pool  boiling  in  horizontal  narrow  annuli  are  reported 
The  effects  of  fluid  properties,  pool  subcooling,  crevice  length,  and  gap  size  on  the 
boiling  behavior  and  the  critical  heat  flux  (CHF)  are  also  studied.  The  CHF  decreases 
with  decreasing  gap  size  or  increasing  length  of  the  annuli.  The  lower  CHF  of  narrow 
crevices  may  be  explained  by  the  thin  film  evaporation. 

.  A  semi-empirical  correlation  is  established  for  the  CHF  of  pool  boiling  in 
horizontal  confined  spaces.  This  correlation  is  compared  with  the  CHF  data  of  the 
present  experiment  Satisfactory  agreement  is  obtained 
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Bond  number  for  the  gap,  defined  in  equation  (2) 

A  constant  defined  in  equation  (1) 

Critical  heat  flux 

Outside  diameter  of  the  heated  tube  [mm] 

Two-phase  friction  factor 

Grashof  number,  g  (Tw  -  T)  d2  /  v[ 

Acceleration  of  gravity  [m/s2 3 
Length  of  the  crevice  [mm] 

Mass  flow  rate,  defined  in  equation  (5)  [kg/s] 

Heat  flux  [kw/m2] 

Critical  heat  flux  [kw/m2 3 

Critical  heat  flux  predicted  by  Kutateladze’s  correlation  for  a  horizontal 
wire,  defined  in  equation  (8)  [kw/m2] 

d,  /  i 

Ambient  temperature  [°C3 

Saturated  temperature  of  the  fluid  [°C3 

T*  ~  Tm  [°C3 

Wall  temperature  of  the  heated  tube  [°C3 
Vapor  velocity  at  the  exit  of  the  crevice  [m/s] 

Width  of  baffle,  defined  in  reference  [7]  [mm] 

Coordinate  measured  in  axial  direction,  defined  in  Figure  11. 


GREEK  SYMBOLS 


Mean  void  fraction,  the  range  between  0  and  1 
Thermal  expansion- coefficient  of  liquid  ['C’l 
Gap  thickness  of  the  crevice  [mm] 

Latent  heat  of  vaporization  [J/kg] 

Viscosity  of  vapor  phase  [kg/ms] 

Viscosity  of  liquid  phase  [kg/ms] 

Kinematic  viscosity  of  liquid  [mVs] 

6f  Angle  measured  from  the  top  of  the  heated  tube  [degree] 

/>,  Density  of  liquid  phase  [kg/mJ] 

Pt  Density  of  vapor  phase  [kg/m*] 

9  Surface  tension  of  liquid  phase  [N/m] 
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INTRODUCTION 


Boiling  at  the  shell  side  of  tube-and-shell  heat  exchangers,  evaporators,  and 
reboilers  in  the  power  and  chemical  processing  devices  received  substantial  attention  in 
the  past  Boiling  may  also  occur  in  the  annular  crevice  between  the  baffle  plate  and 
the  heating  tube.  The  boiling  phenomena  and  the  critical  heat  flux  in  confined  vertical 
annuli  with  closed  bottom  have  been  studied  by  Yao  and  Chang  [1.2].  The  forced 
convective  boiling  in  vertical  confined  annuli  has  been  studied  by  Hung  and  Yao  [3]. 
Presently,  the  boiling  phenomena  and  the  critical  heat  flux  of  confined  annuli  at 
horizontal  orientation  are  studied  in  this  paper  with  specific  application  to  horizontal 
tube-and-shell  heat  exchangers. 

Saturated  pool  nucleate  boiling  on  single  tubes  and  horizontal  plates  has  been 
extensively  studied.  Interest  has  been  extended  to  nucleate  boiling  heat  transfer  in 
restricted  geometries.  Katto  and  Yokoya  et  al.  [4,5.6]  investigated  the  nucleate  pool 
boiling  between  two  restricted  horizontal  plates.  A  movable  artificial  restriction,  which 
is  an  optical  assembly  including  an  interference  plate  and  a  device  for  observation  was 
‘placed  above  the -heated  surface.  It  was  found  that  the  GHF  decreased  with  decreasing 
-gap  size.  At  the  same  heat  flux,  the  wall  temperature  is  invariable  while  she 
interference  plate  is  placed  comparatively  apart  from  the  heated  surface  (  with  gap  size 
larger  than  2.0  mm  ),  but  the  wall  temperature  falls  as  the  gap.  size  is  reduced  and 
thereafter  rises  when  the  gap-  size  is  further  reduced  to  very  .small.  The  dryout 
mechanism  is  suggested  as  a  phenomenon  occurred  when  there  is  an  unbalance  between 
the  consumption  of  the  liquid  by  evaporation  on  heated  surface  and  the  supply  of 
liquid  through  the  intermittent  jetting  of  vapor. 


Improved  heat  transfer  characteristics  with  the  restriction  might  be  attributed  to  an 


increase  in  the  heat  transfer  coefficient  due  to  vaporization  from  the  thin  liquid  film 
on  the  beating  surface  [4,7]  or  increased  bubble  activity  [83.  Jensen,  Cooper,  and 
Bergles  [7]  performed  experiments  of  saturated  water  pool  boiling  at  atmospheric 
pressure  in  horizontal  annuli  utilizing  an  electrically  heated  inner  surface.  Crevice  heat 
transfer  coefficients  were  as  much  as  230  %  greater  than  those  measured  for 
conventional  pool  boiling.  The  increase  in  the  heat  transfer  was  explained  by  the  thin 
film  evaporation.  The  critical  heat  flux  was  found  to  be  directly  proportional  to  the 
gap  size  and  inversely  proportional  to  the  length  of  annulus.  However.  Jensen  et  al. 
[7]  did  not  report  detailed  visual  observations  in  the  narrow  crevice. 

As  for  the  visual  investigation  of  pool  boiling  phenomena,  the  existing  literature  is 
very  limited.  Ishibashi  and  Nishikawa  [9]  studied  the  saturated  pool  boiling  heat 
transfer  in  a  vertical  narrow  annulus  with  both  ends  open.  They  observed  that  there  is 
a  remarkable  difference  of  heat  transfer  between  the  coalesced  bubble  regime  and  the 
isolated  bubble  regime.  These  regimes  were  separated  by  a  critical  gap  size,  which 
varied  with  fluid  properties.  The  isolated  bubble  regime  was  present  in  larger  crevices. 
Below  the  critical  crevice  dimension,  the  coalesced  bubble  regime  was  observed.  Yao 
and  Chang  [1]  presented  a  series  of  systematic  investigations  of  pool  boiling  heat 
transfer  in  vertical  narrow  annuli  with  closed  bottoms.  They  found  that  the  Bond 
number  is  important  in  characterizing  the  boiling  behavior  in  confined  spaces. 
However,  no  systematic  observation  of  pool  boiling  heat  transfer  in  horizontal  narrow 
annuli  can  be  found  in  open  literature. 

The  objectives  of  this  paper  are  to  establish  a  systematic  data  base  of  pool  boiling 
heat  transfer,  in  narrow  horizontal  annuli  with  various  working  fluids  and  provide  the 
explanation  of  the  physical  phenomena:  and  to  develop  a  correction  for  the  CHF 
prediction  in  horizontal  annular  crevices. 


The  liquid  pool  is  established  in  a  Pyrex  tube  of  101.6  mm  inside  diameter  and 
457-5  mm  length  with  both  ends  closed.  The  liquid  level  is  maintained  by  an  equalizer 
connected  to  the  pool,  and  the  temperature  of  the  pool  is  maintained  by  an  immersion 
heater,  which  is  adjusted  by  a  Variac.  A  calibrated  J-type  ungrounded  thermocouple  is 
installed  in  the  pool  to  measure  the  liquid  pool  temperature.  The  system  is  kept  at 
atmospheric  pressure  and  monitored  by  a  pressure  gauge.  The  generated  vapor  in  the 
experiments  is  released  and  condensed  before  flowing  back  to  the  equalizer. 

A  schematic  of  the  test  section  is  shown  in  Figure  L  The  heated  section  is  a 
stainless  steel  304  seamless  tubing  with  0.71  mm  wall  thickness.  25.4  mm  O.D.,  and 
101.6  mm  heated  length.  Both  ends  of  the  heated  tube  are  thermally  shrink-fitted  to  a 
copper  tube  which  has  1.42  mm  wall  thickness.  The  direct  current  passing  through  the 
heated  tube  is  provided  by  two  D.G  power  supplies  of  different  capacities.  The 
current  is  measured  as  the  voltage  drop  with  a  digital  multi-meter  through  a  calibrated 
shunt. 

The  hollow  quartz  cylinders  are  milled  to  an  O.D.  of  63.5  mm  with  the  I.D.  of 
26.04.  27.00,  and  30.56  mm  respectively  to  form  different  gap  sizes  with  respect  to  the 
heated  tube.  Both  the  inside  and  the  outside  surfaces  of  the  quartz  are  polished  to  a 
50-80  finish  to  permit  visual  observation  of  the  annular  crevice.  Four  types  of  annuli 
are  formed  in  the  present  experiments.  Three  of  them  are  76.2  mm  long  with  the  gap 
sizes  of  0.32,  0.80,  and  2.58  mm  respectively.  The  fourth  annulus  is  25.4  mm  long  with 
a  gap  size  of  0.32  mm.  The  geometry  of  the  concentric  annulus  between  the  heated 
tube  and  the  quartz-cylinder  is  maintained  by  two  spacers  located  at  both  ends  of  the 
quartz.  The  horizontal  level  of  the  test  section  is  maintained  during  the  experiment. 
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A  pair  of  movable  thermocouples  are  used  to  measure  the  inside  wall  temperature 
of  the  heated  tube.  Two  J-type  stainless-steel-sheathed,  ungrounded  thermocouples  of 
0.81  mm  diameter  are  pressed  against  the  inner-wall  by  plate  springs  with  a  force  of 
2.6  Newtons.  The  location  of  the  thermocouples  can  be  traversed  axially  or  routed 
circumferentially  at  the  interior  of  the  heated  tube.  The  outer  surface  temperature  of 
the  tube  is  calculated  based  upon  one-dimensional  steady  sute  heat  conduction  [10.11]. 
The  heated  surface  is  polished  before  each  test  with  a  #  320  sandpaper  to  ensure  a 
consistent  condition  of  surface  finish.  The  calibrated  thermocouples  are  connected  to 
an  Accurex  Autodau  Temperature  Logger(  model  Ten/S  )  with  the  accuracy  of  ±0.1°C 

Before  each  experiment  the  test  section  is  preheated  and  preboiled  for  at  least  one 
hour  at  a  low  heat  flux.  The  preheating  for  water  test  required  much  longer  time  than 
other  fluids  to  remove  the  entrapped  air  bubbles.  The  steady  sute  dau  are  uken  after 
3  minutes  for  any  change  of  heat  flux. 

RESULTS  AND  DISCUSSION 

The  experiments  are  performed  with  Freon-113.  Acetone,  and  Distilled  Water  at 
atmospheric  pressure.  The  effects  of  fluid  properties,  gap  size,  crevice  length,  and 
subcooling  on  boiling  curves  and  critical  heat  flux  have  been  investigated. 
Circumferential  wall  temperature  distributions  are  measured.  In  all  the  experiments,  the 
early  dryout  surts  at  the  top-center  point  of  the  heated  tube.  This  is  because  that  at 
very  high  heat  flux  the  top-center  portion  is  difficult  to  be  rewetted  by  incoming 
liquid  [7,12].'  Therefore,  the  critical  heat  flux  in  the  present  experiments  is  measured 
with  the  movable  thermocouple  at  the  top-center  position  of  the  heated  tube. 

In  addition  to  these  measurements,  a  16  mm  color  movie  series  is  edited  as  a 


record  of  the  observations. 


(A).  Characteristics  of  the  Boiling  Curves 


Single-phase  natural  convection.  The  natural  convection  around  a  horizontal 
heated  circular  cylinder  in  the  range  (  104  S  Gr  £  3X10*  )  was  studied  by  Hermann 
[131  and  Jodlbauer  [14].  They  concluded  that  the  Nusselt  number  is  proportional  to 
(Gr)w\  where  Gr  is  based  on  the  diameter  of  cylinder.  That  is.  the  heat  flux  q  is 
proportional  to  (Tw  -  T^)  V4.  Therefore,  the  relationship  becomes 

q  »  C  (  Tw  -  T^  )  *“  (1) 

The  experimental  data  of  single-phase  natural  convection  of  a  single  tube  with 
Freon-113  is  shown  in  Figure  2.  The  C  value  in  equation  (1)  is  found  to  be  0.123  for 
Freon-113. 

The  temperature  distributions  of  Freon-113  at  various  measuring  locations  with  gap 
size  of  0.32  mm  and  CTevice  length  of  76.2  mm  are  shown  in  Figure  3. 

The  location  where  y/L  equals  zero  represents  the  edge  of  the  crevice.  At  the 
same  gap  size  and  0f  (  the  location  of  =  0  is  at  the  top  of  heated  tube  ),  the 
temperature  is  increased  with  the  increasing  y/L,  that  means  the  single-phase  heat 
transfer  is  the  worst  at  the  center  plane  of  the  heated  crevice.  The  temperature  at  the 
bottom  of  heated  tube  is  lower  than  that  at  the  top  of  heated  tube.  Similar  results  are 
observed  for  acetone  and  distilled  water.  It  is  observed  that  the  trend  of  temperature 
distribution  is  almost  the  same  for  various  gap  sizes.  However,  the  smaller  the  gap 
size,  the  higher  the  wall  temperature  for  a  same  heat  flux. 

Boiling  Inception.  When  heat  flux  is  increased  an  overshoot  of  the  wall  superheat 
beyond  saturation  happens  before  the  inception  of  the  boiling.  In  the  present 


experiment  an  interesting  phenomenon  which  should  be  pointed  out  is  the  stability  of 
the  boiling  inception.  Two  different  kinds  cf  boiling  inception  are  observed  :  the  local 
boiling  inception  which  occurs  at  the  top  of  the  heated  tube,  and  the  overall  boiling 
inception  which  occurs  simultaneously  over  the  whole  tube.  Without  observation,  if  the 
temperature  measuring  position  is  at  the  bottom  of  the  heated  tube,  the  local  boiling 
inception  at  the  top  of  the  tube  can  not  be  detected  easily. 

For  the  tests  of  Freon-113  and  acetone  with  2.58  mm  gap  size  and  open  tube  ( 
Le.  single  tube  without  confinements  ),  the  local  boiling  inception  occurs  first,  then  the 
overall  boiling  inception  happens  at  a  higher  heat  flux.  For  the  tests  using  distilled 
water  the  phenomenon  of  local  boiling  inception  is  usually  not  observed.  The  existence 
of  local  boiling  inception  in  horizontal  annular  crevices  could  be  related  to  their  Bond 
number,  which  relates  the  gap  size  and  the  capillary  constant. 


Bo  = 


la/g(PrP)V 


The  range  of  Bond  number  in  the  present  experiments  is  listed  in  Table  1. 


Before  the  boiiing  inception,  the  top  portion  of  the  heated  tube  where  the 
temperature  is  comparatively  high  will  be  more  likely  to  initiate  the  boiling. 


For  the  case  of  Bo 2: 1.0,  after  the  boiling  inception  the  bubbles  are  close  to 
the  spherical  shape  and  they  flow  away  easily.  When  the  heat  flux  is  gradually 
increased,  the  bubbles  may  also  initiate  at  the  bottom  portion  of  the  tube  that  the 
overall  boiling  inception  occurs. 


For  the  cases  of  Bo  <  1.0,  the  generated  bubbles  are  usually  deformed  in  the  gap. 
This  phenomenon  is  also  observed  in  [1].  During  the  boiling  inception  the  bubbles 
expand  into  the  form  of  isolated-dcformed  bubbles  and  propagate  in  the  gap.  Then 
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they  coalesce  over  the  whole  heated  surface.  Therefore,  only  the  overall  boiling 
inception  is  observed. 

Nucleate  pool  boiling.  The  nucleate  boiling  curves  of  the  experiments  are  shown  I 

in  Figures  4.  5,  and  6.  At  a  same  heat  flux,  the  wall  temperature  is  decreased  with 
reducing  gap  sixes  except  for  the  cases  of  acetone  and  distilled  water  with  the  smallest 
gap  of  0.32  mm.  I 

The  increased  heat  transfer  within  the  narrow  crevices  might  be  attributed  to  the 
increased  bubble  activity  [8]  and  the  vaporization  of  the  thin  film  on  the  heated 
surface  [7.91.  In  confined  spaces,  at  a  fixed  heat  flux,  the  masses  of  vapor  generated 
is  constant  so  that  with  decreasing  gap  size  higher  vapor  velocities  are  induced.  With 
the  increased  vapor  velocities,  the  shear  stress  on  the  liquid  film  at  the  heated  surface 
increases  and  the  liquid  film  is  reduced  in  thickness.  Since  the  major  heat  transfer 
resistance  is  the  heat  conduction  across  the  liquid  film,  the  reduced  film  thickness 
increases  the  heat  transfer  coefficient. 

In  the  present  experiments  of  acetone  and  distilled  water  with  gap  size  of  0.32 
mm  and  high  heat  flux,  the  wall  temperature  is  higher  than  that  in  the  cases  of  gap 
sizes  of  0.80  and  2.38  mm  as  shown  in  Figures  5  and  6.  This  is  because  the 
intermittency  of  partial  dryout  and  rewetting  due  to  fluid  oscillation  leads  to  a  wall 
temperature  oscillation  (  the  shaded  area  in  Figures  )  and  a  higher  wall  temperature. 

Similar  phenomenon  has  been  observed  by  Katto  eL  al  [4],  Jensen  et  al  [7],  and  Yao 
et  al  tlj.  . 

*  •  . 

Critical  heat  flux.  At  high  heat  fluxes,  because  of  the  space  restriction,  numerous 
coalesced  vapor  bubbles  occur,  grow,  and  rapidly  spread  over  the  heated  surface.  The 
area  of  heated  surface  is  wet  by  a  thin  liquid  film  while  the  vaporization  takes  place. 

During  the  growth  period  of  the  bubble,  liquid  can  not  be  supplied  to  the  heated 


surface.  The  bubble  gTows  continuously  until  one  of  the  following  two  situations  occurs 
:  the  liquid  film  underneath  the  expanding  bubble  evaporates  completely  and  thereby 
terminates  the  driving  force  for  growth,  or  the  bubble  grows  to  such  a  size  that  it 
extends  beyond  the  edge  of  the  restriction.  Then  the  buoyancy  makes  the  bubble  to 
flow  out  of  the  crevice,  and  the  liquid  flows  into  the  crevice  and  rewets  the  surface. 
After  rewetting,  the  same  cycle  starts  again.  As  the  heat  flux  is  increased,  dry  spots 
begin  to  appear  at  the  end  of  the  cycle.  If  liquid  does  not  rewet  the  surface  within  a 
certain  time,  the  temperature  of  the  dried  portion  may  become  high  enough  to  initiate 
dryout  earlier  in  the  next  cycle.  Therefore,  the  CHF  occurs.  With  the  fluid 
circulation  retarded  when  the  gap  is  reduced,  rewetting  of  the  heated  surface  will  be 
more  difficult  and  the  CHF  becomes  lower.  The  present  CHF  data  at  saturated 
condition  are  listed  in  lable  1. 

(B).  Effect  of  Measuring  Location 


Typical  boiling  curves  measured  at  the  top  and  the  bottom  positions  of  the  heated 
tube  are  compared  in  Figure  7.  This  is  the  case  when  Bond  number  is  less  than  unity 
that  the  overall  boiling  occurs  in  the  gap  at  the  boiling  inception.  Before  the  boiling 
inception  the  heat  transfer  at  the  bottom  of  heated  tube  is  higher.  At  the  overall 
boiling  inception,  the  temperature  decrease  at  the  top  is  more  than  that  at  the  bottom. 
After  the  boiling  inception,  the  heat  transfer  at  the  top  portion  is  dominated  by  thin 
film  evaporation,  and  the  heat  transfer  coefficient  is  therefore  higher  than  that  at  the 
bottom.  As  heat  flux  is  further  increased,  the  thin  liquid  film*  underneath  the  squeezed 
bubbles  may  dry  up  easily.  Finally,  the  permanent  dryout  at  the  top  may  occur  and  an 
early  CHF  is  reached. 


(O.  Effect  of  Crevice  Length 


In  the  experiments  of  Freon-113  with  gap  size  of  0.32  mm.  as  the  length  of 
crevice  is  decreased  from  76.2  mm  to  25.4  ram,  the  overall  boiling  curve  is  changed  as 
shown  in  Figure  8.  The  single-phase  natural  convective  heat  transfer  of  long  crevice  is 
lower  than  that  of  the  short  crevice  due  to  the  higher  resistance  of  the  flow.  After 
the  boiling  inception,  the  local  nucleate  boiling  heat  transfer  at  the  top-center  of  long 
crevice  length  is  higher  mainly  because  the  heat  transfer  at  the  top  is  dominated  by 
thin  film  evaporation.  At  higher  heat  flux,  it  is  difficult  to  rewet  the  top-center 
portion  of  the  long  crevice,  then  the  CHF  of  the  longer  crevice  occurs  at  a  lower  heat 


(D).  Effect  of  Pool  Subcooling 


In  the  subcooled  tests,  the  pool  temperature  is  kept  at  25  °C  subcooling  with  the 
crevice  gap  size  of  0.32  ram.  The  boiling  curve  of  subcooled  condition  is  compared 
with  that  of  saturated  condition  in  Figure  9.  As  expected,  the  subcooling  enhances  the 
boiling  heat  transfer. 


(E).  Boiling  Curve  of  the  Line-Contact  Configuration 


Figure  10  shows  the  boiling  curve  for  the  configuration  of  line-contact  of  the 

f  •  c  •  ■  r  " 

heated  tube  and  the  quartz  shroud  at  the  top  of  the  eccentric  crevice.  When  heat  flux 
is  increased  to  6.05  kw/m2,  a  permanent  dryout  zone  is  observed  at  the  top  contacted 
location.  According  to  the  experimental  results  by  Baum  and  Curlee  [151.  the 
occurrence  of  permanent  dryout  in  some  steam  generators  may  lead  to  local  chemical 
concentration  build-up.  and  the  consequent  corrosion.  As  shown  in  the  Figure  10  the 
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maximum  heat  flux,  which  does  not  lead  to  permanent  dryout  when  a  line-contact 
occurs,  is  much  less  than  the  CHF  of  the  corresponding  situation  with  concentric 
configuration  of  the  annulus. 

(F).  Semi-empirical  CHF  Correlation 


A  correlation  has  been  developed  by  Jensen  et  al.  [7]  relating  the  critical  heat 
flux  to  the  geometry  of  the  horizontal  tube-baffle  crevice  and  fluid  properties.  Using 
a  least-square  fit.  the  following  empirical  equation  was  obtained  : 


The  comparison  of  the  equation  (3)  with  the  present  CHF  data  is  shown  in  Table 
1.  The  prediction  of  equation  (3)  shows  satisfactory  agreement  with  the  present  data 
for  the  small  gap  of  0.32  mm  and  the  short  crevice  length  of  25.4  mm.  but  this 
equation  is  not  adequate  foT  the  long  crevices.  This  is  because  the  test  conditions  in 
[7]  correspond  to  the  cases  of  the  small  gap  size  and  short  crevice  length  in  the 
present  experiments.  In  order  to  predict  the  CHF  in  a  wide  range  of  the  present  data, 
a  semi-empirical  modelling  is  established.  , 


The  schematic  of  the  test  section  is  shown  in  Figure  11-a.  Due  to  the  symmetry, 
the  analysis  has  only  to  apply  for  a  quarter  portion  of  the  whole  annulus.  The  flow 
channel  for  the  present  analysis  is  shown  in  Figure  11-b.  The  gap  size  is  generally 

assumed  much  smaller  than  the  radius  of  heated  tube  such  that  the  curvature  of  the 

•  °  .  * 

crevice  can  be  neglected.  The  flow  channel  of  the  crevice  is  split  into  two  flattened 
regions  as  shown  schematically  in  Figure  11-b  with  the  additional  assumptions  : 

L  The  flow  is  driven  by  the  buoyancy  in  the  Region  I,  and  balanced  by  the 
viscous  drag  force  in  the  Region  II. 


v%  A  _N  .NA, 
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2.  The  fluid  quality  at  the  bottom  of  the  annulus  is  zero.  When  the  critical 
heat  flux  occurs  the  fluid  quality  at  the  upper  exit  of  the  annulus  is  unity. 


3L  The  overall  width  of  the  channel  where  the  vapor  flows  out  at  the  top 
portion  of  the  annulus  (Region  II)  is  d  /2. 


Following  the  assumptions  (1)  and  (2),  the  balance  of  the  buoyant  driving  force 
and  the  frictional  force  can  be  expressed  as 

-  L/2  -  -  2 


-  L/2  p  v 

•v  </>,-/>,>*  di  -  t  I  27  3  £  2  3  (4) 

where  is  the  mean  void  fraction  in  Region  I.  In  this  analysis,  we  assume  =  1/2. 


The  mass  flow  rate  of  the  vapor  at  the  exit  is 


\  (  4  1  1 


e> 


The  critical  heat  flux  can  be  obtained  from  the  conservation  of  energy,  that  is 


•  L  ffd, 

■V  ■  <W ( 5 » <  -T  > 


(6) 


Combining  Equations  (4).  (5).  and  (6),  the  predicted  CHF  qCHF  becomes' 


km 


(7) 


where  qk<(  is  the  critical  heat  flux  predicted  by  Kutateladze’s  correlation  for  a 

t*  t  •  *  •  *  « 

horizontal  wire  in  liquid  pool  [163. 


1/4 


qkB|  =  K,n  \  ptm  t  <r  g  (  p-pt  )  3 

with  the  average  value  of  K,/2  set  at  0.14. 


(8) 


The  F(f)  is  a  function  of  the  two-phase  friction  factor.  The  two-phase  friction 


factor  f  is  dependent  upon  ^  and  p  [17],  and  it  should  be  also  affected  by  Bond 
number  for  narrow  annuli.  Therefore,  the  function  F  can  be  expressed  generally  in  the 
following  non-dimensional  form 


FCf)  =  F{  H  Bo  )  (9) 

The  semi-empirical  function  F  can  be  determined  from  the  best  fit  to  the  present 
data  of  experiments.  That  gives 


F(  .  Bo  )  *  110  (  ^  )  0  3  exp  (  -  0.65BO  )  (10) 

Finally,  the  present  semi-empirical  correlation  is  obtained  from  equations  (7).  (8), 
and  (10). 


The  experimental  data  and  the  prediction  of  this  semi-empirical  formulation  have 
been  compared  in  Table  1.  The  average  error  is  18  %  with  the  maximum  error  less 
than  28  %.  In  the  present  study,  the  maximum  effect  of  the  term  in  the 

correlation  for  various  working  fluids  is  26  %. 


CONCLUSIONS 


Based  on  the  present  study,  the  following  conclusions  may  be  drawn  : 

L  Two  types  of  boiling  inception  are  observed,  i.e.  local  boiling  inception  and 
overall  boiling  inception.  This  phenomenon  can  be  related  to  the  Bond 
number.  Generally,  no  local  boiling  inception  is  observed  when  the  Bond 
number  is  less  than  unity. 

2.  After  the  overall  boiling  inception,  the  boiling  heat  transfer  at  the  top 
portion  of  the  heated  tube  is  possibly  dominated  by  thin  film  evaporation. 
During  all  the  experiments,  the  early  critical  heat  flux  occurs  at  the  top- 
center  point  of  the  heated  tube. 

3.  Sustained  flow  oscillation  is  observed  in  the  tests  of  0.32  mm  gap  size. 
Occasionally,  the  wall  temperature  of  this  narrow  crevices  at  nucleate  boiling 
is  higher  than  that  of  larger  gaps.  This  is  due  to  the  intermittency  of  partial 
dryout  and  rewetting  which  leads  to  an  obvious  temperature  oscillation  and 
the  higher  wall  temperature  in  narrow  crevices. 

4.  Boiling  phenomena  of  an  eccentric  annulus  with  the  configuration  of  line- 
contact  at  the  top  of  heated  tube  are  investigated.  A  permanent  dryout  zone 
occurs  at  the  top  portion  of  the  heated  tube  at  low  heat  flux. 

5.  The  CHF  decreases  with  decreasing  gap  size  or  increasing  length  of  the 
annuli.  The  mechanism  of  CHF  in  confined  spaces  is  related  with  the  dryout 
of  the  thin  liquid  film  on  the  heated  surface. 
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Table  1  :  The  Range  of  Bond  Number  and  the  CHP  Data  in  the  Present  Experiments 


Tube 


Operating  Pressure 
Pool  Temperature 


O.D.  :  25.4  mm 
Thickness  :  0.71  mm 

Materia]  :  Stainless  Steel  304  Seamless  Tubing 
1  atm 

25  °C  Subcooling  -  Saturation 


Woring  Fluid 

L 

S 

Bo 

^CHF 

Pred. 

^CHF 

%  Error 

(mm) 

(mm) 

(kw/m2) 

(kw/m2) 

of  [II] 

[11 

[in 

Freon-113 

76.2 

0.32 

0.30 

112 

17.5 

8.8 

-27.8 

FTeon-113 

76.2 

0.80 

0.74 

27.8 

36.6 

26.1 

-  6.3 

Freon-113 

76.2 

158 

138 

63.7 

98.2 

51.9 

-18.6 

Freon-113 

25.4 

0.32 

0J0 

35.7 

41.5 

45.7 

28.4 

Acetone 

76.2 

0.32 

0.20 

16.1 

30.8 

14.3 

-11.3 

Acetone 

76.2 

0.80 

0.50 

55.0 

64.6 

46.5 

‘  -15.4 

Acetone 

76.2 

2.58 

1.60 

110.5  . 

173.3 

131.7 

19.2 

Acetone 

25.4 

0.32 

0.30 

94.7 

73.2 

74.2 

-21.6 

Distilled  Water 

•  •  .  * 

76.2 

0.32 

0.13 

33.9 

58.4 

4L1 

21.4 

113  :  CHF  predicted  by  the  equation  (3).  Ref.  [73. 

CHI  :  CHF  predicted  by  equations  of  (7),(8),and  (10). 
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Chapter  2 

Heat  Convection  and  Fluid  Flow 
in  Tube  Bundles  at  Very  Low  or 
Very  High  Reynolds  Numbers 
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C.  D 


VC1 


Nu  .  Nu 

h  b 

Nu,.  Nu, 


cylinder  radius 


Momenclature 


,*A,a 


fluid  envelope  radius.  ( - ) 

V 

coefficients  in  eq.  (11) 
coefficients  in  eq.  (11) 
constants  in  eq.  (34) 


nondimensional  coefficient.  Re  Pr  St 

Dh 

cylinder  diameter,  2a 
hydraulic  diameter,  defined  in  eq.  (2) 
hydraulic  diameteT,  defined  in  eq.  (3) 
Fanning  fraction,  defined  in  eq.  (21) 


heat  transfer  coefficient 


thermal  conductivity  of  fluid 
local  and  average  Nusselt  number 
local  and  average  Nusselt  numbers 


pressure 


Peclet  number,  Pr  Re„ 

d 

Prandtl  number,  v/a 


radial  coordinate 


Reynolds  number,  Ud/v 

Reynolds  number,  UD„/v 

modified  Reynolds  number,  defined  in  eq,  (24) 

solidity.  *dJ/[4XTXL] 

Stanton  number 


nondimensional  local  fluid  temperature, 
local  fluid  temperature 


T  -  T 


T  -  T 

b  * 


-.v -v- .v.v.v.v  ••.v 


■r.  + 


wall  temperature 

fluid  bulk  mean  temperature 

free  stream  velocity 

flow  velocity  at  minimum  cross-section  in  bundles 

radial  component  of  fluid  velocity 

6  component  of  fluid  velocity 

longitudinal  pitch 

transverse  pitch 

nondimensional  variable, 

thermal  diffusivity  of  fluid 


ratio  of  cell  radius  to  tube  radius,  — 

a 


voidage,  1  -  S 

nondimensional  stream  function, 
coordinate 
function  of  solidity. 


UaX 

f  "  1 


(l+yff4)  In  f  +  (l-/?4) 


viscosity  of  fluid 
kinematic  viscosity  of  fluid 
variable,  defined  in  eq.  (30) 
radial  shear  stress  component 
tangential  stress  component 


stream  functiton 
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ABSTRACT 


The  empirical  formulation  of  Biery  tl],  which  employed  geometric  transformation, 

for  the  heat  transfer  in  tube  bundles  at  high  Reynolds  number  and  cross  flow  has  been 

—  .6 

extended.  A  simple  heat  transfer  correlation  can  be  deduced  as  Nu  =  0.366  Re, 
Pr,/3.  This  is  because  the  homogenous  turbulent  mixing  in  the  bundles  at  high  Reynolds 
number  at  cross  flow  the  local  heat  transfer  dominates  that  detail  geometry  of  bundles 
is  no  more  important 


At  very  low  Reynolds  number,  the  flow  streamlines  in  a  bundle  are  symmestric 
with  respect  to  a  plane  separating  the  upstream  and  downstream  side  of  a  tube.  The 
flow  field  can  be  approximated  by  the  "free  surface  model"  togather  with  creeping  flow 
approach.  The  calculations  show  that  the  pressure  drop  is  dependent  upon  on  the 
solidity  (1/fl1)  and  transverse  pitch  to  diameter  ratio  (XT/d)  of  a  tube  bundle.  The 

O 

Fanning  friction  factor  is  described  as  f  =  — — (ft  /  H  )0  M  ,  where  Re*  =  Re  [  In  p1 
1  -  Re* 

+  -  HI  -  d/XJ*1.  This  derived  equation  can  predict  experimental  data  of 

1  +  fir 

friction  factor  within  ±  20%  for  Reynolds  number  Re,  less  than  200.  For  the 

9 

condition  of  high  Peclet  numbeT  the  average  Nusselt  number  on  a  tube  is  derived  as 

—  A4  ~  1 

Nu  =  1.185  X,/3  Re1/3  Prl/3  (u./ u  )0M  ,  where  X  = - .  Compared 

4  4  (l+A4)  in  A2  *  (l-A4) 

with  the  existing  heat  transfer  data  in  the  range  of  Red  <  200,  a  slightly  modified 
coefficient  is  proposed  to  fit  the  data  in  ±  30%.  This  gives  Nu,  =  0.83  X1'3  Re,,/3  Pr1'3 


e 
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INTRODUCTION 


Extensive  studies  [1-10]  has  been  conducted  for  the  heat  transfer  in  tube  bundles 
at  cross  flow  due  to  its  wide  range  of  applications.  Large  quantity  of  data  was 
presented  in  reference  [10]  for  heat  transfer  and  pressure  drop  over  smooth  tube 
bundles  at  high  Reynolds  numbers.  From  the  designer’s  view  point,  there  is  a  need  of 
a  general  correlation  on  the  prediction  of  heat  transfer  coefficients  and  friction  factor 
of  cross  flow  through  tube  bundles.  Recent,  Biery  [1]  proposed  a  new  method  by 
transforming  the  actual  tube  bundle  to  an  equivalent  equilateral  triangular  pitch  infinite 
smooth  tube  bank  (ETP-I-STB)  based  on  the  same  hydraulic  diameter  to  tube  diameter 
ratio.  The  heat  transfer  results  is  presented  in  the  form  of 


Nu_  =  C  Re  0  6  Pr1 

Dh  h  Dh 


It’s  found  that  the  transformation  had  the  property  of  bringing  all  of  the  various 
configurations  of  tube  bundles  onto  a  single  curve  of  Ch  vs.  (XT/d)^  (transverse  pitch 
to  diameter  ratio  of  ETP-I-STB).  The  results  were  compared  with  published 
experimental  data  for  bundles  with  pitch  to  diameter  ratios  range  from  1.25  to  2.50  and 
Reynolds  numbers  Re0  above  4000.  The  prediction  of  Cb  for  staggered  array  tube 

ll 

bundles  appeared  to  be  quite  satisfactory  but  an  error  of  20  percent  too  high  existsfor 
in-line  tube  bundles.  The  hydraulic  diameter  in  Biery’s  study  [1]  was  defined  based  on 
the  minimum  flow  cross  section  area  in  tube  bundle  as  that  of  Kays  &  London  [9]. 


4  (  minimum  flow  cross  section  area  )  (  longitudinal  pitch  ) 
(  heat  transfer  area  ) 


While  in  some  of  the  early  literatures  [2-93  the  hydraulic  diameter  was  defined  based 
on  the  volume  of  fluid  in  tube  bundles. 


:• 


\ 


(3) 


n  =  4  (free  volume  ) 

*  (  wetted  perimeter  ) 

At  low  Reynolds  number  condition.  Gunter  and  Shaw  [5]  recommended  an 
empirical  correlation  for  the  pressure  drop  of  cross  flow  through  tube  bundles.  The 
single  correlation  was  generally  in  good  agreement  for  many  tube  layouts  of  interest, 
but  was  considerably  lower  in  prediction  for  the  square  staggered  tube  arrangements 
[3,5].  On  the  other  hand.  Ishihara  and  Bell  Ill]  solved  the  Navier-Stokes  equation, 
numerically  for  in-line  square  tube  bundles  with  pitch  to  tube  diameter  ratios  range 
from  1.20  to  125  and  Reynolds  number  from  1  to  100.  Their  results  for  1.25  and  1.50 
pitch  to  diameter  ratios  are  in  well  agreement  with  experimental  data.  But  no  general 
correlation  of  friction  factor  was  proposed. 

At  the  low  Reynolds  number  condition,  the  heat  transfer  data  of  staggered  tube 
bundles  have  been  conducted  by  Whitaker  [12].  The  correlation  is  well  satisfied  with 
the  data  of  most  staggered  bundles  but  can  not  predict  the  data  of  in-line  bundles  well. 

At  high  Reynolds  numbers,  it  is  likely  that  the  tubulent  mixing  is  so  substantial 
that  even  a  simpler  formulation  than  that  of  Biery  [1]  can  be  generated.  On  ther 
limit  of  very  low  Reynolds  numbers,  the  creeping  flow  type  approximation  can  be  used 
to  derive  the  formulation  of  friction  factor  and  heat  transfer.  This  is  the  objective  of 
this  paper  to  present  the  results  of  these  studies  at  extreme  Reynolds  numbers. 

HIGH  REYNOLDS  NUMBER  CONDITIONS 

In  the  transformation  approach  of  Biery  [1],  the  coefficient  C  in 

b 


bundles  was  correlated  as  : 


IV.W. 


V  •-  \  V  *. 


Ch  =  0.2818  (XT/d)A  -  0.1282  - 


0.08263  [  1.2  -  (XT/d)Al 
[(XT/d)A  -  0.8924  ] 


1  uh 

Where  (  XT/d  =  -y(3.628  -  +  l)l/2  is  the  transverse  pitch  to  diameter  ratio  of  the 

equivalent  equilateral  triangular  pitch  infinite  smooth  tune  bundle  (  ETP-I-STB).  the 
Reynolds  number  Re_  is  based  on  the  hydralic  diameter  D  of  eq.  (2)  and  the  flow 

Dh  h 

velocity  at  minimum  flow  cross  section  in  bundles. 


In  fact,  the  eq.  (1)  and  eq.  (4)  can  be  rearranged  into  a  simple  expression  as 


Nu.  =  C  Re  0  4  Pr1/J  ( - )°  4 

4  -  «*  Dh 


=  f  [  0.09533  «•  0.1409  (  3.628  —  1  ),a] 

t  A 


0.034308 

[0.3628  D /d  +  1]  1/2  -  0.9848 

o 


)  Re  °‘  Pr1,J  [ — ]04 


Where  Nud  is  the  average  Nusselt  number  around  a  tube  in  the  bundles. 
Examination  of  eq.  (5)  shows  an  interesting  result  that  C  (d/D  )0  4  is  nearly  a 

a  h 

single  value  of  0.366  for  the  tube  spacing  range  from  1.25  to  2.50  as  shown  in  Fig.  1. 
Thus  eq.  (5)  can  be  simplified  as  : 


Nu  =  0.366  Re.0*  Pr,/J 

o  0 


Hence,  equation  (6)  gave  a  simple  expression  on  the  heat  transfer  coefficient  of  tube 
bundles  in  spite  of  the  tube  spacing  and  tube  layout  in  this  described  range  of 
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conditions.  Eq.  (6)  is  compared  with  various  experimental  data  in  Fig.  2.  It  can  be 
seen  that  the  data  tend  to  fall  within  i  20%  of  eq.  (6)  for  Re^  greater  than  4000.  It's 
believed  that  at  high  Reynolds  number  conditions  the  violent  tubulent  mixing  exists  in 
the  bundles  such  that  the  local  effect  at  tube  wall  dominate  and  the  detail  geometry  of 
bundle  is  no  more  important.  However  this  phenomena  breaks  down  at  intermediate 
Reynolds  number  conditions  as  shown  in  Fig.  2  especially  for  in-line  tube  bundles. 

LOW  REYNOLDS  NUMBER  CONDITIONS 


When  the  Reynolds  number  Rej  is  low.flow  field  can  be  well  characterized  by 
creeping  flow  approximation.  At  this  conditions  the  flow  streamlines  are  symmetric 
with  respect  to  a  plane  which  passes  through  the  tube  and  normal  to  the  incoming 
stream. 

Various  investigations  [13,14]  of  creeping  flow  through  packed  beds  and  fluidized 
beds  has  been  conducted  by  using  "free  surface  model".  The  application  of  this  model 
reduces  the  complex  and  difficult  problem  of  the  multiparticle  assemblage  by  a  set  of 
concentric  cells.  The  model  assumes  that  in  this  cell  a  spherical  particle  is  enclosed  by 
a  hypothetical  fluid  envelope.  The  envelope  has  zero  shear  stress  at  its  outer  surface. 
Each  cell  contain  the  same  amount  of  fluid  as  the  relative  volume  of  fluid  to  particle 
volume  in  the  entire  assemblage.  The  entry  and  exit  effect  are  neglected  in  this  model. 
The  fluid  mechanics  study  of  creeping  flow  in  tube  bundles  using  the  "free  surface 
model"  has  been  performed  in  reference  [16].  However  no  correlation  on  pressure 
drop  were  recommended  and  no  analysis  on  heat  transfer  was  performanced.  In  fact, 
the  objective  of  this  study  is  to  provide  this  information  to  the  readers. 


ANALYSIS 


Fig.  3  shows  the  diagrams  of  the  typical  staggered  array  and  in-line  array  tube 
bundle  with  the  flow  normal  to  the  tubes.  In  the  following  section  the  fluid  flow  and 
heat  convection  will  be  analyzed  with  the  following  assumptions. 


L  The  flow  is  steady  state  with  constant  fluid  properties. 

1  The  fluid  is  at  laminar,  creeping  flow  with  negligible  energy  disspation. 

3.  The  bundles  are  infinite  with  negligible  entry  and  exit  effects.  The  tube  end 
effect  are  also  negligible. 

4.  Cell  model  is  representative  for  all  the  geometry  of  bundles.  This  assumption 
will  be  justified  by  the  comparision  of  the  present  results  to  reported 
experimental  data. 


With  the  use  of  the  "free  surface  model",  the  hydrodynamic  problem  of  the 
bundles  is  reduced  to  one  of  a  cylinder  with  radius  a  and  surrounded  by  a  cylindrical 
envelope  of  fluid  with  a  free  surface  at  radius  b  as  shown  in  Fig  4.  The  fluid 
envelope  moves  as  cross  flow  at  a  constant  velocity  U,  which  is  the  velocity  at  the 
mininum  cross  section  in  the  tube  bundle.  The  ratio  of  the  cell  radius  to  the  tube 
radius  .  /},  is  related  to  the  voidage  «  of  the  bundle  by 


fi  =  ( - - - )  1/2  (7) 

1  -  t 


1 

The  >S  is  related  to  the  solidity  of  a  bundle  S  as  S  =  — .  The  voidage  of  tube 

bundle  can  be  also  related  to  the  pitch  to  diameter  ratio  as 


JLX,  -  »rd2  /  4 

-LI -  (8) 


v. 


*.  N  ,  «  _  * 
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Therefore  one  gets  : 


This  is  the  relationship  between  the  tube  bundles  pitch  XT,  XL  with  the  unit  cell 
radius  b  in  the  free  surface  model. 


FLUID  FLOW 


Formulation.  For  simplicity  of  derivation,  the  cylinder  of  radius  a  is  assumed  to 
moves  perpendicular  to  its  axis  in  a  fluid  cell  of  radius  b  as  shown  in  Fig.  4.  The 
momentumn  equation  of  creeping  flow  with  inertia  terms  neglected  can  be  written  in 
term  of  stream*  function. 


V4  f  =  0 


The  general  soluton  for  eq.  (10)  is  found  as  : 


f’(r  ,9)  -  sin#  [  Cr3  +  Drlnr  +  Er  +  —  ] 


The  boundary  conditions  in  the  cell  model  are  [12]: 


V  *  U 


V  =  U  COS# 
r 


=  U  sin  5 


at  r  =  a 


Theses  conditions  are  based  on  the  no-slip  conditions  on  the  cylinder  together  with 
the  zero  shear  stress  conditions  at  the  outer  surface  of  the  unit  cell.  They  provides  the 


simultaneous  equations  to  evaluate  the  constants  in  equation  (11). 


a*  U 

C  - - 

2  (a4+b4)  ln(a/b)  +  0.5  (a4  -  b4) 


D  * 


-  U 

ln(b/a)  ♦  a4/(a4+b4)  -  (1/2) 


U  tab 

ta(b/a)  +  (a4-b4)/2(a4+b4) 


F 


-  U  a2  b4 

2(a4+b4)  ln(b/a)  ♦  (a4  -  b4) 


(13) 


Therefore,  the  stream  function  is 


X  r  U  sin0  r2  f? a2 

=  - - -  (  —  ♦  2  (1+f)  ln(b/r) - —  )  (14) 

(/-  1)  X  a2  r2  1 


To  obtain  the  velocity  component  Vf  and  on  the  basis  that  the  inner  cylinder 
is  stationary  and  that  the  fluid  envelope  is  moving  with  the  constant  velocity  U,  one 
simply  needs  to  add  -U  cos#  to  and  U  sin#  to  which  can  be  derived  from  eq. 
(14).  Therefore,  the  correspondingly  stream  function  becomes 

f  -  f'  -  r  U  sin# 

X  U  r  sin#  r2  /a2 

- -  -  (  +  2(1+/)  ln(a/r)  -  (  1-/) - —  \  (15) 

(/-l)  X  a2  r2  1 


and 


V  a 


X  U  cos# 


-yy-  (  +  2(1+/)  ln(b/r)  -  )  -  U  cos# 


-  X  U  sin# 


V#  = 


j  5intr  ? 

— —  (  3  -r-  -  2(1+/)  On  r  +1)  +  2(1+/) 


(/- 1  ) 


lnb 


■&T-  )  +  U  sin# 


(16) 


/  -  i 

Where  X  =  - ,  which  decreased  with  increasing  solidity  1/#*. 

(1+/)  In/  +  (1-/) 


Drag  Force.  Attention  is  drawn  here  upon  the  magnitude  of  the  pressure  gradient 
in  tube  bundles  which  is  of  great  importance  in  practice.  The  method  given  in  [153 


was  used  to  calculate  the  drag  force  on  a  cylinder  from  the  knowledge  of  stream 
function  by  integrating  the  surface  force  around  the  surface. 


Fd  =  £  (<r  cosfl  -  cr^sinS)  a2  dd 


(17) 


where  <rf  and  <t  q  represent  the  radial  and  tangential  stress  component  and  defined 
as  : 


,  «v, 

a  =  -  p  +  2u - 

r  3r 


i  av 


r  3r 


L) 


(18) 


Substituting  V  ,  from  eq.  (16)  into  eq.  (18)  and  (17).  eventually  one  gets 


•  8  f  ^  a  U  (^  1)  1 
(/  -  1) 


(19) 


The  force  per  unit  length  associated  with  a  single  cell  may  be  related  to  the 
pressure  gradient  in  tube  bundles.  Hence, 


y*  _  dP 

wb2  dx 


(20) 


The  Fanning  friction  factor  f  is  defined  in  term  of  the  pressure  gradient  and 


hydraulic  diameter  [15]  by  the  following  equation. 


(-  -jf-)  -  *  r  .  — 


As  mention  early,  two  different  definition  of  hydraulics  diameter  were  made  by 
several  investigators.  Use  the  definition  of  eq.  (2)  following  Kays  and  London  [10]. 


4  X,  \ 

7  <-r-»>-r 


(#dJ/4) 


In  the  present  study,  the  definition  of  hydrodynamic  diameter  in  eq.  (22)  is  found 
to  be  superior  than  the  one  following  eq.  (3)  in  correlating  the  data.  The  equation 
(19)  to  (22)  can  be  rearrange  to  give 


Where  the  modified  Reynoplds  number  Re*,  which  is  related  to  the  bundle 
geometry,  is  defined  as  : 

(1-/94) 

Re*  =  Re,  [  In  f  * -  3  <  1  -  d/X  )  * 1  (24) 

d  (1-A4)  ^ 


To  account  for  the  property  change  of  fluid  at  heating  or  cooling,  a  correction 
factor  [23  /  — —  \  0,4  is  applied  in  eq.  (23).  Finally,  the  friction  factor  is  becomes 


V. 


8  (/+!)  (1  -  d/X^ 


Z.  r  (  —)  014 

Re,  [  2(yS4+l)  ln^  ♦  il-f)  3  '  • 


Where  //b  and  are  the  fluid  viscosity  evaluated  at  Tfe  and  respectively. 


Discussion.  The  calculated  Fanning  friction  factor,  eq.  (25),  is  plotted  in  Fig.  5 
versus  tube  Reynolds  number  for  flow  across  in-line  square  tube  bundles  at  different 
pitch  to  diameter  ratio.  Comparision  is  made  with  Bell’s  numerical  solution  [10]  and 
published  experimental  data  [3,6,7].  It  appears  that  the  Fanning  friction  factor 
increased  as  pitch  to  diameter  ratio  decreased  and  decreased  with  increasing  Reynolds 
number.  Present  results  are  in  well  agreement  with  the  experimental  data  and  Bell’s 
solution  at  p/d  -  L25,  1.5  at  various  Reynolds  numbers.  It’s  also  noticed  that  the 
creeping  flow  results  is  only  slightly  effected  by  inertia  for  Re,  below  200. 

Further  comparision  of  the  pressure  drop  data  for  bundles  with  various  tube 
arrangement  is  shown  in  Fig.  6  where  the  calculated  fricton  factor  is  plotted  versus  the 
modified  Reynolds  number  Re*.  All  the  available  data  tend  to  lie  within  ±  20%  of  the 
present  study  for  Re*  <  400.  When  the  modified  Reynolds  number  is  greater  than  400, 
turbulence  probably  appears  in  various  degree  in  the  tube  bundles,  and  the  creeping 
flow  approach  under-predicts  the  pressure  drop  in  tube  bundles.  Although  eq.  (25)  is 
more  complicated  than  Gunter  and  Shaw’s  correlation  [5],  the  present  correlation 
appears  to  be  a  somewhat  better  one  because  it  fits  all  the  data  for  Re,  <  200,  while 
Gunter  and  Shaw’s  correlation  predicted  50%  low  values  for  staggered  bundles  with 
square  tube  layout  [9]. 


HEAT  TRANSFER 


V.V.V.V/.'; 


*irViinrw 
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3lT  1  AT 

Du*  to  the  thin  boundary  layer  (r-a)/a  «  1,  - »  —  The  second  term  at 

dr2  r  ®r 

the  right  side  of  eq.  (27)  is  negligible.  The  steady  state  convective-diffusion  equation, 
therefore,  can  be  rewritten  in  the  form  of 


at  _ 

Hd~ 


4  (2  sin*)1'1 

Pe4  X 


at 


*  (  17  ) 


(29) 


where 


2  U  a 


is  the  Peclet  number 


Using  the  new  variable 
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eq.  (29)  becomes 
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Eq.  (31)  can  be  further  simplified  by  introducing  the  similarity  variable  Z. 
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(2  sin*)1'2  y  (Pe„  X)1/J 

_  Q 

a  (4  (2)m)m 
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The  governing  eq.  can  be  trasnsf erred  into  the  ordinary  differential  equation. 


1L  ♦  ±  7}  JL  -  ft 
dzl  3  Z  dZ  0 


Eq.  (33)  can  be  integrated  twice  : 


Tz  dZ  +  c. 


The  constant  c2>  c}  can  be  evaluated  from  boundary  conditions 


(0  at  r  =  a,  t  *  0  .  therefore 

Cj  *  0 

(ii)  at  r  =  b.  (  9  ->  oo.  Z  ->  oo)  .  t  =1  ,  therefore 


The  second  boundary  condition  is  rewritten  as  that  the  stream  function  at  the 
outer  surface  of  the  unit  cell  approach  infinite.  This  is  reasonable  since  the  quantity 
(b-a)  is  always  much  greater  than  the  thermal  boundary  layer  thickness  at  High  Peclet 
number  conditions.  Finally,  one  can  get 


\  -  Tw  (9/4)  H4/3) 


Nusselt  Number.  The  local  Nusselt  number  based  on  the  hvdraulic  diameter  is 


obtained  from  the  equaton. 
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Substituting  eq.  (35)  into  eq.  (36)  yields 
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Therefore,  the  local  Nusselt  number  based  on  the  tube  diameter  can  be  written  as 


Nud  (i 9 .  Ped)  = 


(16  Pe  X),/3  (sin*)1'2 

_  O 
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In  this  equation  the  dependence  of  the  heat  transfer  coefficient  on  the  hydraulic 
diameter  of  cross  flow  in  tube  bundles  disappears.  The  average  Nusselt  number  Nu  is 

d 

evaluated  by  integrating  the  local  Nusselt  number  around  the  tube  . 


Nu.  =  1.185  Re : 1/3  Xl/3  Pr1/3 
•  « 


The  correction  fractor  (  —  \  0,14  is  also  proposed  here  to  account  for  the 

'  /i*  ' 

properities  change.  Therefore,  the  average  Nusselt  number  becomes  : 


0.14 


(39) 


Nu4  =  1.185  Re4,/J  Pr,/3  (  —  ) 

'  M  * 


Equation  (39).  therefore,  give  an  analytical  relationship  between  the  heat  transfer 
coefficient  and  Reynolds  number  for  tube  bundles  of  any  regular  geometry  with  the 
solidity  (  )  of  the  tube  bundle  as  a  parameter. 

P 


Discussion.  The  present  solution  is  restricted  to  low  Reynolds  number  and  high 

NUj 

Peclet  number.  Fig.  7  shows  the  comparison  of  eq.  (39)  and  data  in  term  of  - 

pr'/3  hJ 

versus  X  Re,  for  various  tube  arrangement  We  may  also  note  that  the  data  of  which 

Q 

natural  convection  is  significant  [16,17]  are  excluded  in  present  analysis.  The  eq.  (39) 
over-predicts  the  data,  this  is  possible  because  the  range  of  Reynolds  number  of  the 
data  base  is  larger  than  the  range  of  the  creeping  flow  that  the  flow  streams  are,  in 
fact  non-symmtric  for  the  upstream  and  downstream  portion  with  respect  to  the  tube. 

The  analysis  is  assuming  creeping  flow  gives  a  symmetric  stream  lines  which,  therefore, 
over-predicts  the  heat  transfer.  In  order  to  fit  the  data  base  in  this  Reynolds  number 
range,  it  is  suggested  that  the  constant  in  eq.  (39)  is  reduced.  That  suggests 


Nu4  =  0.83  Re*  1,1  Pr,,J  X1'3  (  —  )0M  (40) 


All  the  data  fall  within  ±  30%  of  this  semi-empirical  correlation.  Eq.  (40) 

contains  a  term  X  which  is  dependent  upon  the  tube  spacing.  The  increasing  of  the 

1 

tube  spacing  results  in  lower  heat  transfer  coefficient  In  other  words,  the  solidity  ( — 

P 

)  plays  an  important  role  in  the  performance  of  cross  flow  in  tube  bundles  at  low 
Reynolds  number  high  Peclet  number  conditions. 
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To  compare  in  a  comon  base  it  is  found  that  Whitaker's  correlation  [12]  can  be 
reformed  in  a  similar  form. 


Nud  =  1.526  Red,/J  -  1)“2/3  Pr1/J  (  —  )° 14 

P  w 


Equation  (41)  and  (40)  are  compared  with  experimental  data  in  Fig.  8  for  the 
effects  of  /S2.  Both  eqs.  show  the  tendence  of  increasing  the  heat  transfer  coefficient 
by  increasing  the  solidity  (1/y?2)  of  the  bundles.  The  comparison  indicates  that  the 
present  semi-empirical  correlation,  eq.  (40),  is  slightly  better  than  the  empirical 
correlation  of  Whitaker's. 


CONCLUSION 


The  empirical  formulation  of  Biery  [1],  which  employed  geometric  transformation, 
for  the  heat  transfer  in  tube  bundles  at  high  Reynolds  number  and  cross  flow  has  been 
extended.  A  simple  heat  transfer  correlation  can  be  deduced  as  shown  in  eq.  (6).  This 
is  because  the  homogenous  turbulent  mixing  in  the  bundles  at  high  Reynolds  number  at 
cross  flow  the  local  heat  transfer  dominates  that  detail  geometry  of  bundles  is  no  more 
important  However  this  phenomena  breaks  down  at  intermediate  Reynolds  number  as 
shown  in  Fig.  2.  Eq.  (6)  give  a  satisfied  prediction  of  ±  20%  error  compared  with  the 
experimental  data  for  Reynolds  number  above  4000. 


At  very  low  Reynolds  number,  the  flow  streamlines  in  a  bundle  are  symmestric 
with  respect  to  a  plane  separating  the  upstream  and  downstream  side  of  a  tube.  The 
flow  field  can  be  approximated  by  the  "free  surface  model"  togather  with  creeping  flow 
approach.  The  calculations  show  that  the  pressure  drop  is  dependent  upon  on  the 
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INTRODUCTION 


Boiling  heat  transfer  in  tube  bundles  has  been  extensive  investigated  for  many 
years  due  to  various  engineering  application.  This  problem  is  important  in  the  design 
of  nuclear  reactor,  evaporator,  condenser,  chemical  reboiler  and  many  major  components 
in  chemical  and  thermal  power  plants. 

Boiling  heat  transfer  can  be  divided  into  two  categories  :  pool  boiling  (  where 
the  tube  bundle  is  submerged  in  an  stagnant  pool  )  and  flow  (forced  convection  ) 
boiling  (  in  which  fluid  flow  in  the  heated  tube  bundle  ).  Estimation  of  boiling  heat 
transfer,  CHF  and  pressure  drop  on  the  outside  of  tube  bundles  is  one  of  the 
important  factor  in  successful  design  of  heat  exchanger  or  steam  generator. 

Most  of  previous  investigation  [1-63  however,  were  concerned  with  the  axial  flow 

along  tube  bundles.  As  a  matter  of  fact,  in  many  thermal  system  the  steam  generator 

and  evaporator  contain  tube  bundles  which  are  oriented  horizontally.  The  flow 
direction  is  normal  to  the  axial  of  the  tubes  and  boiling  occurs  in  the  shell  side  as 

shown  in  Fig.  1.  The  tubes  can  be  arranged  in  in-line  or  staggered  array  with  a 

particular  pitch  to  diameter  ratio,  which  determines  the  tube  to  tube  distance.  Steam 
and  water  flows  by  the  gravity  induced  natural  convection.  However,  at  any  position  in 
the  bundle  the  boiling  heat  transfer  on  the  tube  is  mainly  dependent  upon  the  local 
flow  velocity  and  the  local  quality.  The  study  of  heat  transfer  in  a  particular  region 
of  a  large  bundle  can,  therefore,  be  performed  in  a  smaller  bundle  at  a  forced 
convective  boiling  condition  with  the  same  velocity  and  quality  environment 

Reliable  prediction  of  heat  transfer,  CHF  and  pressure  drop  is  vital  in  optimiming 
the  heat  exchanger  design.  Unfortunately,  the  knowledge  of  the  hydraulic  and  boiling 


heat  transfer  in  tube  bundles  at  cross  flow  condition  is  very  inadequate.  Extensive 
reviews  on  heat  transfer  of  single-phase  cross  flow  in  tube  bundles  are  presented  in 
reference  [7-16]  In  two-phase  flow  case,  the  vapor  bubble  distribution  affects  the  flow 
field  and  it  can’t  be  handled  by  the  well  known  knowledge  of  single-phase  flow  in 
tube  bundles. 

Fundamental  understanding  and  parametric  effect  has  been  fairly  well  established 
for  the  boiling  heat  transfer  and  CHF  of  a  single  tube  at  cross  flow  condition  in 
reference  [17-21].  The  major  difference  between  the  single  tube  heat  transfer  and  the 
beat  transfer  on  a  tube  in  a  heated  tube  bundle  is  the  different  flow  field  and  quality 
environment  near  the  tube  surface. 

The  available  information  of  cross  flow  boiling  heat  transfer  in  tube  bundles  is 
mainly  oriented  from  the  experience  of  the  reboiler  design  in  chemical  process  [21-26] 
and  the  evaporator  design  [27-29].  No  adequate  data  exist  in  the  literature  for  the 
effect  of  quality  and  velocity  on  the  critical  heat  flux  in  tube  bundles.  The  peripheral 
variation  of  heat  transfer  coefficient  on  the  heated  tube  in  single  tube  or  bundles  are 
also  not  available  in  published  literature. 

It's  also  interesting  to  study  the  staggered  array  bundle  because  the  flow  streaming 
phenomena  (  see  Fig.  l.b  )  may  not  occurs  in  staggered  array.  Therefore,  the  heat 
transfer  coefficient  could  be  higher  because  most  of  the  tube  surface  area  would  be 
active.  The  comparision  of  boiling  heat  transfer  in  staggered  array  bundle  with  that  in 
in-line  array  bundle  is  still  not  clear  up  to  now.  For  a  better  understanding  of  boiling 
heat  transfer  in  tube  bundles,  a  systematical  study  is  required  to  bridge  the  well  known 
single  tube  knowledge  to  tube  bundles.  Twro  tubes  or  three  tubes  study  [28-29]  has 

I 

served  as  a  base  for  the  study  of  tube  bundles. 


The  objective  of  the  present  investigation  is  to  bridge  the  knowledge  of  heat 
transfer  on  single  tube  to  the  tube  bundles,  and  to  compare  the  heat  transfer  of  tube 
bundles  with  staggered  array  and  in-line  array.  The  effect  of  the  thermal  environment 
of  a  tube  in  a  bundle  can  also  be  revealed  by  comparing  the  heat  transfer  results  of 
only  one  tube  heated  vs.  all  tubes  heated,  and  single  phase  flow  vs.  two  phase  flow. 
With  the  better  understanding  on  the  boiling  heat  transfer  and  two  phase  flow,  the 
computation  of  the  thermal  hydraulics  in  steam  generator  bundles  may  become  possible. 


The  major  effort  in  the  past  year  is  to  fabricate  and  assembly  the  test  section  and 
the  test  loop.  At  the  same  time,  the  analysis  of  single  phase  foced  convection  in 
various  tube  bundles  at  very  low  or  very  high  Reynolds  number  has  been  completed. 
This  analysis  was  presented  in  the  chapter  2  of  this  report  The  experimental  set  up 
will  be  reported  in  the  following  section. 


EXPERIMENTAL  APPARATUS 


•  The  major  components  of  the  experimental  apparatus  are  :  a  Freon-113  loop,  test 

section,  power  supply  system,  instrumentation  system  and  photographic  equipment  All 
of  these  components  will  be  described  in  detail  in  the  following  sections. 

Freon  Loop 


The  Freon-113  loop  in  Thermal  Science  Lab.  at  Carnegie-Mellon  University  was 
•  modified  to  achieve  a  wide  range  of  operation  condition.  The  schematic  of  the 

modified  loop  is  presented  in  Fig.  2  with  its  specifications  described  in  Table  1. 


All  parts  of  the  main  flow  system  are  constructed  of  type  304  stainless  steel  and 


insulated  with  fiber  glass.  Freon-113  is  circulated  by  a  Crane  stainless  leakage-proof 
dynapump,  model  JB-3K,  which  provides  a  discharge  pressure  of  80  psig  at  30  GPM. 
After  the  pump,  the  flow  is  divided  into  the  by-pass  line  and  the  test  section  line.  In 
the  test  section  line  ,  fluid  flows  through  the  turbine  flow  meter,  the  preheaters,  the 
main  regulating  valve,  and  the  test  section.  After  the  test  section,  the  fluid  merges 
with  the  flow  from  the  by-pass  line.  The  mixed  flow  then  goes  through  the  condenser, 
degassing  tank,  and  then  flows  back  to  the  pump.  Two  Grainger  filters  aTe  installed  in 
parallel  at  by-pass  suction  line  of  the  pump. 

The  flow  rate  through  the  test  section  is  adjusted  by  the  valve  at  by-pass  line  and 
the  valve  at  the  inlet  of  the  test  section.  The  desired  system  pressure  is'  maintained  by 
the  Nitrogen  gas  in  the  accumulator.  Pump  operation  is  controlled  also  by  the 
regulating  valve  at  the  outlet  of  the  pump.  The  main  regulating  valve  (  1.5  inches 
diameter  ),  which  is  located  at  the  entrance  of  the  test  section,  is  used  to  regulate  the 
two  phase  quality  into  the  test  section.  Therefore,  the  inlet  condition  to  test  section 
will  be  controlled  by  varying  the  power  delivered  to  Freon-113  in  the  preheaters  and 
the  throttling  at  the  main  regulating  valve  before  test  section. 

The  preheater  system  consists  of  four  Chromalox  heaters  which  has  the  capability 
of  producing  a  total  of  30  KW  of  power.  The  first  and  second  preheater  are  Model 
ARMTO-3605T2,  each  with  a  rated  output  of  6  KW  at  240  V.  The  third  and  fourth 
preheater  are  Model  MTO-390A  of  9  KW  each.  Each  preheater  is  able  to  provide  a 
continuously  variable  power  using  variacs.  The  preheaters  are  arranged  in  a  way  to 
prevent  the  burn  out  by  trapped  bubbles  on  the  heaters  The  second  preheater  is  laid  at 
the  top  of  first  preheater.  The  fourth  preheater  is  installed  above  the  third  one  with  a 
10  degree  inclination.  A  separate  valve  on  top  of  the  second  preheater  is  also  used  to 
bleed  the  trapped  air  in  pt  eaters. 
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The  condenser  is  of  a  shell-and-tube  type  with  Freon-113  flowing  in  the  shell 
side  and  cooling  water  at  the  tube  side.  The  valve  on  the  water  line  is  able  to  control 
the  temperature  of  Freon-113  before  the  pump  inlet  After  the  condenser.  Freon-113 
enters  the  degassing  tank,  which  is  located  at  the  highest  point  of  system. 

A  photograph  of  the  loop  assembly  is  shown  in  Fig.  3.  The  range  of  parameters 
for  operation  conditions  in  present  study  is  shown  in  Table  2. 


Test  Section 


The  test  section  under  investigation  contains  the  bottom  plenum,  test  bundle  and 
top  plenum.  The  tubes  in  the  test  bundle  can  be  reused,  the  test  section  casing,  inlet 
plenum  and  exit  plenum  are  the  same  in  all  the  tests. 


The  schematics  of  the  test  section  is  indicated  in  Fig.  4  with  its  photograph 
presented  in  Fig.  5.  The  typical  cross-section  of  the  test  section  is  shown  in  Fig.  6. 


Three  different  bundles  are  fabricated  for  geometries  A,  B  and  C  as  shown  in  Fig. 
7.  The  geometries  are  also  specified  in  Table  3.  Geometry  A  is  mainly  used  to  study 
the  heat  transfer  and  CHF  of  a  single  tube  in  a  line  of  tubes.  The  heated  tube  or  the 
unheateds  tube  can  be  drawn  out  of  the  test  section  from  the  back  plate  to  produce 
different  kind  of  spacing.  Geometry  B  is  a  in-line  tube  bundle  with  pitch  to  diameter 
ratio  of  1.5.  Geometry  C  is  a  staggered  tube  bundle  with  1.5  pitch  to  diameter  ratio 
in  both  transversal  and  longitudinal  directions.  -  Theses  two  bundles  are  design  to  give 
the  same  solidity,  which  means  the  same  number  of  tubes  per  unit  cross  section  area  of 
the  bundle,  in  order  to  compare  their  heat  transfer  behavior  at  a  same  base. 


As  described  in  Table  3,  the  test  bundle  consists  of  an  upstream  unheated  zone,  a 
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heated  zone,  and  an  exit  unheated  zone.  The  tubes  in  the  unheated  zone  are  solid 
Aluminum  tubes,  while  the  tubes  in  the  heated  zone  are  304  stainless  steel  seamless 
tubes.  The  inlet  unheated  zone  is  used  to  stabilize  the  flow  stream  and  to  allow  a 
smoth  flow  transition  from  the  inlet  plenum  to  test  bundles.  It  is  also  designed  to 
compare  the  pressure  drop  with  that  in  the  heated  zone. 

The  casing  of  the  test  bundle  is  made  of  rectangular  Aluminum  duct  with  660.4 
mm  length.  The  cross  section  of  flow  duct  is  61.6  mm  by  85.725  mm.  The  duct  is 
vertically  orientated  in  the  flow  loop  with  Freon-113  flowing  against  gravity.  Glass 
windows  (  3/4  inches  thick  pyrex  No.  7740  glass  )  are  installed  in  several  locations  to 
enable  visual  observation,  lighting  and  photographic  recording  of  flow  pattern.  The 
front  window  openings  are  installed  with  the  2.5  inches  diameter  pyrex  glasses.  The 
side  window  openings  are  installed  with  1.75  inches  diameter  pyrex  glasses.  A  1.5  inch 
level  shift  between  two  sides  windows  is  made  to  obtain  a  better  lighting  condition. 
The  1/16  inches  thick  neoprene  rubber  gaskets  between  the  glasses  and  the  Aluminum 
casing,  and  the  Teflon  tape  are  used  to  avoid  the  direct  contact  of  pyrex  glass  and 
Aluminum  casing  and  to  seal  the  test  section. 

The  changeable  stainless  stell  back  plate  (3/4  inches  thick  )  is  attachecd  to  the 
Aluminum  casing.  The  back  plate  is  fitted  with  electrically  heated  tubes  of  19.05  mm  ( 
3/4  inches  )  diameter,  and  0.508  mm  (  0.020  inches  )  wall  thickness.  The  detail  of  the 
non-instrumented  heating  tube  is  illustrated  in  Fig.  8.  The  electrically  heated  tubes  are 
silver-soldered  onto  two  copper  ends.  These  copper  pieces  and  the  central  copper  rod  ( 
0.281  inches  diameter  )  serve  as  the  electrical  current  passage.  The  central  copper  rod 
is  silver-soldered  on  the  copper  head  at  one  end  and  screwed  to  the  copper  buss  bar  at 
the  open  end.  The  Teflon  disk  mounted  at  the  end  of  the  copper  pieces  is  used  as  the 
electrical  insulation. 


Each  test  bundle  contains  an  instrumented  tube  in  the  heated  zone.  The 


peripheral  temperature  on  this  tube  is  measured  by  two  J-type  ungrounded  stainless 
steel  sheathed  thermocouples  of  1.588  mm  (  1/16  inches  )  diameter.  In  the 
instrumented  heating  tube,  the  Micarta  tube  which  guard  the  central  copper  rod  has  a 
slot  (  1/8  inches  wide  )  cut  through  it  to  provide  access  for  thermocouples  and  ceramic 
insulator.  Fig.  9  shows  the  detail  of  the  instrumented  heating  tube  assembly.  Pressure 
taps  and  thermocouples  are  also  installed  at  several  locations  in  the  test  chamber  to 
provide  the  flow  field  information. 


Both  the  inlet  plenum  and  the  exit  plenum  are  made  of  stainless  304  pipe  of  101.6 
mm  (  4  inches  )  diameter,  50.8  mm  (  2  inches  )  length.  The  inlet  plenum  contains  a 
perforated  plate  to  homogenize  the  inlet  vapor  and  liquid  distribution  of  the  test 
bundle.  The  exit  plenum  contains  flow  straightener  to  prevent  the  uneven  flow 
distribution  near  the  top  of  the  bundle. 


Power  Supply 


The  D.C  electrical  power  for  the  test  heating  tubes  is  provided  from  a  small 
welder  and  four  large  Westinghouse  power  supply  in  the  Thermal  Science  Laboratory. 
The  small  welder  supplies  the  finely  controlled  power  up  to  40  V,  300  Amp,  while  the 
large  Westinghouse  power  supply  gives  up  to  100  V  and  1500  Amp. 


The  D.C  power  from  the  Westinghouse  welders  is  transported  to  the  power  pannel 
through  two  2000  copper  welder  cables  and  then  to  the  copper  buss  bar  of  the  heating 
tube  in  the  test  bundle.  A  current  shunt  is  installed  in  series  with  the  test  section  such 
that  the  high  current  input  to  the  test  tube  bundle  can  be  accurately  measured.  For 
CHF  experiment,  two  current  shunts  are  used  to  measure  the  different  heat  fluxes 
supplied  to  the  instrumented  tube  and  surrounding  tubes. 


Instrumentation 


All  the  measured  temperatures  in  the  test  section  are  recored  by  the  Accurex 
Autodata  logger  model  Ten/5.  The  Autodata  Ten/5  is  a  microprocessor-controlled  data 
acquisition  system  which  can  scan  30  channels.  It  converts  the  output  from 
thermocouples  to  precisely  scaled  DC  voltage  for  measurement  and  displays  the  data  in 
engineering  units.  A  limited  mathematics  package  that  permits  computation  in  real 
time.  The  alarms  are  also  provided  as  a  warning  system  in  CHF  studies.  The  data 
processing  is  programmed  with  a  connected  Zenith-19  cathode  ray  terminal.  The 
photograph  of  the  instrumentation  rack  is  shown  in  Fig.  10. 

The  system  temperatures  in  the  loop  are  measured  by  J-type  (  1/16  inches  ) 
thermocouples  and  read  directly  from  Autodata  Ten/5  and  digital  temperature  indicator 
(  Omega  Model  403A-C  )  through  the  selector  switches.  Two  digital  multimeters  ( 
Model  HP-3465A  )  are  also  used  in  conjunction  with  two  compensators  (  Omega  Model 
MCJ-J  )  to  obtain  more  accurate  temperature  measurements.  Two  moveable 
thermocouples  as  shown  in  Fig.  9  are  used  to  measure  the  inside  wall  temperature  of 
the  heated  tube.  The  junction  of  the  J-type  stainless  steel  sheathed  unground 
thermocouples  of  0.794  mm  diameter  are  pressed  against  the  inner  wall  by  a  plate- 
spring.  ■* 

The  power  to  the  heating  tubes  is  obtained  by  measuring  the  current  through  the 
heating  tubes  or  the  voltage  drop  across  the  precision  stnnL 


The  flow  rate  through  the  test  section  is  measured  with  a  turbine  flow  meter  ( 
Model  80-0750,  manufactured  by  Engineering  Measurement  Company  Inc  ).  A  standard 
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preamplifier  is  supplied  with  the  meter.  The  range  of  the  output  from  the  amplifier  is 
4  to  40  mAmp  corresponding  to  1.5  to  34  GPM.  The  current  output  from  the  meter  is 
converted  to  mVolt  (  8  to  80  mv  )  and  is  fed  to  a  Leeds  &  Northrop  Model 
Speedomax  H  chart  record.  The  line  response  of  the  meter  can  be  described  as  : 


Q  (GPM)  =  0.468  V  (mv)  -  3.75 


The  system  pressures  are  measured  from  0.32  cm  diameter  taps  located  at  several 
locations  in  the  loop  and  in  the  test  section.  The  pressure  are  recorded  with  two 
bourdon  pressure  gages  0  -  500  psi  in  0.5  psi  and  1.0  psi  sub-divisions.  The  pressure 
drop  through  the  tube  bundles  is  also  obtained  by  a  calibrated  pressure  transducer  ( 
Model  PB  415D-13  )  which  is  connected  to  a  chart  recorder  (  Sanborn  Model  321  ). 
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Table  1  Specification  of  Freon-113  Loop 


Item  No. 

Description 

1 

Main  regulating  valve 

2 

Test  section 

3 

Ball  valve 

4 

Condenser 

5 

Global  control  valve 

6 

Degassing  tank 

7 

Filter 

8 

Pump 

9 

Turbine  flow  meter 

10. 

Preheater  1 

11 

Preheater  2 

12 

Preheater  3 

13 

Preheater  4 

14 

Freon-113  supply  tank 

15 


Freon-113  accumulator 


Table  2  Range  of  Experimental  Conditions 


Working  fluid 

Heating  tube 

O.D.  of  heating  tube 

Wall  thickness  of 
heating  tube 

Test  section 

Channel  flow  area 

Heat  flux 

Inlet  fluid  quality 

Superficial  liquid  velocity 

Volume  flow  rate 

Mass  flow  rate 


Freon  -  113 
seamless  304  SS  Tube 

19.05  mm 

0.508  mm 

61.6  x  85.725  x  660.4  mm 
52.82  cm2 
0  -  537  KW/m2 
0  -  0.3 

0  -  0.34  m/sec 
0-30  GPM 
0  -  490  Kg/m2 
0-40  psig 


System  pressure 


Triangular  Array 
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